Oblique dark solitons in supersonic flow of a Bose-Einstein condensate 

G.A. E110 A. Gammal^lil and A.M. Kamchatnov^ 

^ Department of Mathematical Sciences, Loughborough University, Loughborough LEll 3TU, UK 

^ Instituto de Fisica, Universidade de Sao Paulo, 05315-970, C. P. 66318 Sao Paulo, Brazil 

^Institute of Spectroscopy, Russian Academy of Sciences, Troitsk, Moscow Region, 142190, Russia 

(Dated: February 8, 2008) 

In the framework of the Gross-Pitaevskii mean field approach it is shown that the supersonic flow 
of a Bose-Einstein condensate can support a new type of pattern — an oblique dark soliton. The 
corresponding exact solution of the Gross-Pitaevskii equation is obtained. It is demonstrated by 
numerical simulations that oblique solitons can be generated by an obstacle inserted into the flow. 
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Introduction. It is known that the nonlinear and dis- 
persive properties of a Bose-Einstein condensate (BEG) 
can lead to the formation of various nonlinear structures 
(see, e.g., [l[). Until recently, most research has been 
focused on experimentally observed vortices and bright 
and dark solitons. Furthermore, the formation of dis- 
persive shock waves in BEGs with repulsive interactions 
between atoms was considered theoretically in [2, l3| and 
studied experimentally in rotating [4] and non-rotating 
[5| condensate, where it was shown that dispersive shocks 
are generated as a result of the evolution of large distur- 
bances in the BEG. However, another important type 
of nonlinear structure, namely a spatial dark soliton, can 
also be realized in a BEG. The first experimental evidence 
of their generation has recently appeared [6|. In fact, the 
existence of oblique spatial solitons in a BEG has a nat- 
ural physical basis if the Gherenkov generation of disper- 
sive sound waves by a small obstacle in the supersonic 
flow of a BEG is considered and the effect of increasing 
the size of the obstacle (i.e. the amplitude of the waves) 
is determined. Evidently, along with dispersion, nonlin- 
ear effects become equally important at finite distances 
from the obstacle, so that the Gherenkov cone breaks-up 
into a spatial structure consisting of one or several dark 
solitons. Such a structure represents the dispersive ana- 
log of the well-known steady spatial shock generated in 
the supersonic flow of a viscous compressible fluid past 
an obstacle. In this sense, it is the spatial counterpart 
of the one-dimensional expanding dispersive shock ^j- 
[5| generated in the evolution of large disturbances in a 
BEG. In the simplest case, the nonlinear wave structure 
would consist of a single spatial dark soliton given by a 
steady solution of the equations governing the BEG flow. 
Motivated by this physical consideration and the results 
of experiments [61], in this Letter we shall develop the 
theory of spatial dark solitons in the framework of the 
Gross-Pitaevskii (GP) mean field approach. 

Basic equations. The dynamics of a BEG is described 
to a good approximation by the GP equation pj] 



where ?A(r) is the condensate order parameter and g is an 
effective coupling constant, with g = ATrh^as/m, Os be- 
ing the s-wave scattering length and m the atomic mass. 
Here ^(r) denotes the potential of the external forces 
acting on the condensate, for example the confining po- 
tential of the trap and/or the potential arising due to 
the presence of an obstacle inside the BEG. When the 
"obstacle" is formed by a laser beam and the flow occurs 
due to the free two-dimensional expansion of the BEG 
the trap potential should be set equal to zero for the free 
expansion of a BEG, and far enough from the obstacle we 
can neglect the obstacle potential as well. Also, we are 
interested in steady flows, that is, we suppose that the 
parameters of the flow change on a timescale much slower 
than the transient timescale for the establishment of the 
wave pattern of interest. To this end, we seek solutions 
of Eq. ID) with T^(r) = of the form 



il){v) = \/n{v) exp [ - 
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where n(r) is the density of atoms in the BEG, u(r) 
denotes its velocity field and /x is the chemical poten- 
tial. It is now convenient to introduce the dimensionless 
variables f = r/V2^, h = n/no, u — u/c^, where no 
is a characteristic density of atoms, equal to their den- 
sity at infinity, ^ = hj \J2mnQg is the healing length and 
Cs = h/^/2m£^ is the sound velocity in a BEG of density 
no. Substituting Eq. ([2]) into ^ and separating real and 
imaginary parts we obtain a system of equations for the 
density n{x, y) and the two components of the velocity 
field u= {u{x,y),v{x,y)), 
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where we have omitted tildes for convenience. If we 
l^) restrict our consideration to the vortices- free potential 



flows with vanishing curl of the velocity field, 



0, 



(4) 



then the second and third equations in ([3]) can be in- 
tegrated once to give a generalization of the Bernoulli 
theorem to dispersive 2D hydrodynamics: 



-(u^ + v^)+n + ^ (nl + ^1)-^ {n^x + nyy) 



= const. 



(5) 

Equations ([4]), ([5|) and the first equation of ^ comprise 
the system governing the BEC potential flow. 

Our aim now is to flnd the solution of this system under 
the conditions that the BEC flow is uniform at infinity 



n = 1, u = M, V — at 



oo, 
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where M denotes the ratio of the asymptotic velocity of 
the flow to the sound velocity, i.e. the Mach number. 

Oblique dark soliton solution. Let us seek a solution of 
the form n = n{6), u — u{9), v = v{9), where 9 = x — ay 
and a denotes the slope of the soliton center location 
with the y axis. Substitution of this ansatz into ^ and 
the first equation of ([3]) , followed by a simple integration, 
yields expressions for the velocity components in terms 
of the density 



M{l + a^n) 



aM{l - n) 



(7) 



where the integration constants were chosen according to 
condition ©. Then substitution of ([7]) into (O, with a 
proper choice of the constant in the right-hand side, leads 
to the equation 

i(l + a^){n1 - 2nngg) + 2n^ - {2 + p)n'^ +P = 0, (8) 

where 

p^M^il + a^). (9) 

It is easily checked that Eq. ([5]) has the integral 

l{l + a')nl^{l^nf{n~p) (10) 

where, again, the integration constant is chosen in ac- 
cordance with condition ^. Simple integration of this 
equation finally yields the desired solution in the form of 
a dark soliton for the density 



ni9) = 1 
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(11) 



The velocity components can then be found by substitu- 
tion of this solution into Eqs. ([7]). The inverse half-width 
of the soliton in the x-direction is 



1-p 



(12) 



Thus formulae (|lip and (|12p give the exact dark spa- 
tial soliton solution of the GP equation. We shall call 
it "oblique" because it is always inclined with respect to 
the direction of the supersonic flow. Numerical solutions 
below show that such solutions are stable for M > 1. 

Small amplitude Korteweg-de Vries (KdV) limit. As 
is clear from dill) the small amplitude limit is achieved 
when 1 — p ^ 1. Then the parameters a and p can be 
expressed in terms of k and M from Eqs. ^ and (1121) as 
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The density proflle in this limit becomes then familiar 
KdV soliton 



n = 1 — 
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where k < l/M. 

Note that the slope a = \J M"^ — 1 corresponds exactly 
to the Cherenkov cone, so that in its vicinity the small 
amplitude solitons are located inside it. This approxima- 
tion corresponds to the KdV limit of the potential-free 
GP equations ([3]). Indeed, if we assume the series expan- 
sions (e ^ 1) 



n — \ + erii 



u = M + eui 



V = evi 



and introduce the scaling of the independent variables 
C — e^/^{x — ay), t = e'^^'^y, then standard reductive 
perturbation theory leads to the KdV equation for the 
density disturbance ni 



drfli 
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M^ 
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with the well-known soliton solution of this equation 
equivalent to dH]), p^ . 

Nonlinear Schrodinger (NLS) equation limit. Another 
important limit corresponds to large slopes a^ S> 1. For 
this limit we introduce the parameter A as p = M"^ jc? — 
A'^, that is a = ±M/\. The soliton solution ([TT|) can then 
be approximated as 



1-A2 



(17) 



cosh^[^/T^^{y ± \x/M)] ■ 
This is exactly the solution n = |$p of the NLS equation 



i^T + '^YY -2|*r* = 



(18) 



for the complex variable ^ = ^/?i exp I i J v{Y' , t)dY' j , 

where T ~ x/2M and Y = y. This equation was derived 
in jTj from the GP equations ([3]) for the highly supersonic 
(M » 1) flow of a BEC past a slender body. 

Generation of oblique solitons in a BEC. Let us now 
consider the supersonic flow of a BEC past an obstacle. 



If the obstacle is small (e.g. an impurity), then linear 
sound waves are generated at finite distances which form 
a Cherenkov cone Q. Large obstacles generate spatial 
dispersive shocks which can be viewed as trains of inter- 
acting dark spatial solitons inside the Cherenkov cone. 
The theory of the generation of spatial dispersive shocks 
has been developed in much detail for supersonic flows 
past a slender body when such a flow can be described 
by the KdV equation [9| . Analogous theory for the NLS 
equation case was developed in [7|. However, in real ex- 
periments the obstacles cannot be considered as slender 
bodies and the flow is not highly supersonic, hence fully 
nonlinear solutions of the GP equation, such as Eq. (fTTj) . 
should be used for the quantitative description of spatial 
dispersive shocks in a BEC. Here we shall use numer- 
ical solutions of the GP equation to demonstrate that 
the structures generated by an obstacle inserted into the 
supersonic BEC flow indeed contain the oblique dark soli- 
tons given by Eq. pTj) . 

To make this process clearer, numerical solutions of 
the time-dependent GP equation ([1]), expressed in non- 
dimensional variables as 



«V't = -T^ii'xx + i^yy) + Vix,y)^p + IV;] 



(19) 



were studied, where V' — i^/y/no and t = {gno/h)t, with 
the other variables defined as in Eq. ([2]). From now on, 
the tildes will be omitted. The potential V{x,y) corre- 
sponds to the interaction of the condensate with the ob- 
stacle. Since the detailed behavior of the potential should 
not be critical for the formation of solitons far from the 
obstacle, it can be safely modeled by an impenetrable 
disk. In our simulations its radius was set to r = 1 in our 
non-dimensional units. Initially at < = there is no dis- 
turbance in the condensate, so that it is described by the 
plane wave function tpix,y)\^^Q = exp(zMa;) correspond- 
ing to a uniform condensate flow. To be specific, let us 
take M = 5. Several stages of the numerically calculated 
BEC density evolution are shown in Fig. 1. 

It can be clearly seen how a pair of oblique solitons 
is gradually formed behind the obstacle. Their length 
grows with time and, except in the vicinity of the end 
points, the density distributions do not demonstrate any 
vorticity, which agrees with our assumption of potential 
flow Q. However, the flow cannot be considered as sta- 
tionary and potential near the end points. This is mani- 
fested by the presence of vortices behind the end points 
of the spatial solitons. One may interpret these vortices 
as a "vortex street" [lOJ radiated by spatial dark solitons. 
However far enough away from these vortex end points, 
the flow can be considered stationary. The establish- 
ment time of the stationary profile can be estimated by 
the soliton width divided by the sound velocity, which 
is t ~ 2 in our solutions. This is much less than the 
value i = 20 for the last plot of Fig. 1. The parameter p 
from Eq. © was calculated using the value of the slope 




FIG. 1: Emergence of a pair of oblique dark solitons after 
"switching on" supersonic flow {M = 5) past a disk-shaped 
impenetrable obstacle of radius r = 1 located at {x — 0,y = 
0). The direction of the flow is from left to right. Density 
plots are shown for three times t = 5, i = 10 and t — 20. The 
dark structures correspond to oblique dark solitons, which in 
turn generate the vortex streets near the end points. 



a inferred from the numerical solution. A comparison of 
the theoretical profile of the oblique dark soliton given 
by Eq. (fTT|) with the corresponding part of the density 
profile in the full numerical solution is shown in Fig. 2. 
The excellent agreement between these two profiles con- 
firms that the line patterns in Fig. 1 are indeed oblique 
dark solitons generated by the obstacle. This agreement 
also justifies the assumptions made in the derivation of 
the analytic solution (fTTj) . Note that, along with the 
soliton, a small amplitude dispersive wave packet is gen- 
erated, which corresponds to the "non-solitonic" part of 
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FIG. 2: Cross sections of the density distributions for x — 
20 (dashed hne), a; = 60 (sohd line) and j/ > obtained 
from numerical solution of the GP equation (|19|l . These are 
compared with soliton profiles (lllf) with slope a = 10 shown 
as functions of y at the same values of a; (a; = 20 corresponds 
to "crosses" and a; = 60 to "circles"). 



the density perturbation induced by the obstacle. This 
wave packet spreads out as distance from the obstacle 
increases and eventually fades away. 

In the above simulations, the parameters of the ob- 
stacle have been chosen so that only a single soliton is 
generated at each side of the obstacle. However, with 
the increase of the size r of the obstacle the number of 
solitons is also expected to increase. This effect is demon- 
strated in Fig. 3 in which two symmetric fans of solitons 
can be seen behind the obstacle. As expected, the depth 
of the dark solitons grows with the increase of the slope a 
with respect to the y axis. Thus, the oblique dark solitons 
can be viewed as "building blocks" in more complicated 
patterns arising in supersonic flows of a BEC. This fig- 
ure also demonstrates the robustness of the phenomenon 
for different obstacle parameters. The whole structure 
in Fig. 3 represents a pair of dispersive shocks generated 
in the supersonic flow of a BEC past an obstacle. Such 
shocks were considered in [7| in the limiting case of highly 
supersonic flow {M :^ 1) and slender obstacles (a :^ 1). 
Our present numerical simulations show that spatial dis- 
persive shocks represent a general phenomenon caused by 
the interplay of dispersive and nonlinear effects described 
by the full GP equation. We speculate that it can be also 
be observed in flows of a BEC past 3D obstacles. 

The above theory is based on the supposition that the 
flow of a BEC is homogeneous and uniform enough for 
periods of time sufficient for generation of oblique dark 
solitons. Let us indicate here the physical conditions 
which should be satisfied for obtaining such a flow in 
experiment. To be definite, we consider the example of 
2D flow in "pancake" geometry, that is the condensate 
is supposed to be confined in axial direction by a strong 
harmonic potential. Let the transversal frequency ujj_ 
of the potential be small enough so that the (radial) 



FIG. 3: Density plot at the moment t = 30 for the supersonic 
flow (M — 5) past a disk-shaped impenetrable obstacle with 
radius r = 5 located at {x = 0,y = 0). 



Thomas-Fermi density distribution is applicable. We 
consider expansion of the BEC after switching off the 
potential (see, e.g. [11[). The asymptotic as i :^ l/i^± 
solution of the cylindric GP equation for r <ti luj±t where 
I is the radius of the BEC before its release from the 
trap, is given by n(r,t) = constant/(wj_i)^, u{r,t) = r/t, 
that is the density is practically uniform but varying 
with time. In the same approximation, the local heal- 
ing length ^ and the local sound velocity Cg are given by 
^ — ht/ml, Cs = l/{\/2t). Hence, the local Mach number 
M = ^/2r/l. Thus, we get a supersonic flow past obsta- 
cle if we place it at the distance d > l/V^ from the axis. 
Now, the flow can be considered as uniform if the size 
of the obstacle (scaled as healing length ^ for a chosen 
moment of observation) satisfles the condition £,/d <C 1. 
For d ^ I this gives t ^ mP/H, i.e the flow past an ob- 
stacle is asymptotically uniform for wj^ ^ t <C mP /h. 
At last, the characteristic time of establishing the soliton 
distribution, ^/cs ~ m^'^ /h, obviously satisfles the above 
inequality since ^ ^ Z, so that the flow can be considered 
as quasi-stationary. At the same time, our numerical sim- 
ulations show that oblique solitons are generated even in 
non-uniform and nonhomogeneous flows past obstacles, 
that is the phenomenon is very robust with respect to 
change of parameters of the flow. 

To summarize, we have found an exact oblique dark 
soliton solution of the GP equation and demonstrated 
numerically that such solitons can be generated by ob- 
stacles inserted into the supersonic flow of a BEC. 
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